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We show that Gordon's optical metric on a curved spacetime can be generahzed to include absorp- 
tion by allowing the metric to become complex. We demonstrate its use in the realm of geometrical 
optics by giving three simple examples. We use one of these examples to compute corrected distance- 
redshift relations for Friedman-Lemaitre-Robertson- Walker models in which the cosmic fluid has an 
", associated complex index of refraction that represents grey extinction. We then flt this corrected 

^sL ■ Hubble curve to the type la supernovae data and provide a possible explanation (other than dark 

' energy) of the deviation of these observations from dark matter predictions. 
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I. INTRODUCTION 

' Gordon [ij made the interesting observation that any solution to Maxweh's equations in a curved spacetime filled 
\^ with a fluid whose electromagnetic properties can be described by a real permittivity e{x) and a real permeability 
I fi{x), or hence a real refraction index n{x) = ^/e/^i can be found by solving a slightly modified version of Maxwell's 
rS |. equations in a related optical spacetime with vacuum values for the permittivity and permeability, i.e., with e{x) = 1 
' and /i(x) = 1, see Eq. 0. That is to say there is a one to one relation between solutions of Maxwell's equations 
^ , in these two spacetimes, the physical spacetime with an index of refraction and its physical metric, and the optical 
' spacetime with n — I and its optical metric. For traveling electromagnetic waves, the modified geometry of the 
. optical metric accounts for the decrease of the wave speed which in physical spacetime is actually caused by an index 
of refraction n > 1. Until now the open question was whether or not absorption could be incorporated into the 
optical metric. In this paper we show that if the Maxwell field is a monochromatic wave the optical metric can be 
^ . modified to account for absorption as well as refraction. That is, a solution to a slightly modified set of Maxwell's 
— ' equations, without refraction or absorption, in the optical spacetime gives the appropriately absorbed and refracted 
t — , wave in physical spacetime. This conclusion extends to the superposition of multi-frequency waves as long as the 

' optical properties are frequency independent. 
^""^ , We assume that at some given frequency range the fluid's electromagnetic properties are linear and isotropic 
Cf^ ■ relative to the fluid's unit 4-velocity Uau'^ = —1, and can be summarized by a complex permittivity e{x°') and a 
' complex permeability /i(a;") defined on the four dimensional spacetime manifold. Following Q we write a complex 
0^ , refraction index N(x'^) as 

q. ^ ' 

j> ! N ^ ^/e^= n + iK, (1) 

' where n and k are respectively the real and imaginary parts. The real electric field/magnetic induction bivector ^ 
Fab{E, B) = Re{F^l^} that represents a traveling monochromatic wave has a geometrical optics expansion (see 0) of 
■ ■ ■ ' the form 

pXo ^ g,:S/Ao (^^^ ^ ^^^^ ^ OiXo')) , (2) 



and satisfies the homogeneous Maxwell equation^ 



dlaFt^^O, (3) 
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in both the physical and optical spacetimes. We designate the constant expansion parameter of geometrical optics by 
Xo because we use its value to adjust the wavelength. It's positioning as a superscript or subscript is for convenience 
only. The real part of S{x°') is the usual eikonal function which determines the surfaces of constant phase for the wave. 
The Aab term represents the usual amplitude of the geometrical optics approximation and the Bab term is its first 
order correction (covariant components Fab, Aab, and Bab are identical in the physical and optical spacetimes). The 
constitutive relations for the contravariant components of the real displacement/magnetic field bivector in physical 
spacetime H°'''{D, H) = Re{/7^^} are given by 

Ht-\Ft^\-r-rF^s, (4) 

where the optical metric gab of Gordon [l[ (which now becomes complex) has been used to raise the covariant indices 
of F^^^ to produce F^\ The complex optical metric is related to the real physical metric gab = —UaUt + g±ab by 

gab = (1 )UaUb + gab = -^UaUb + g_Lab, (5) 

with inverse 

5'''' =: (1 - e/i)M°u'' + g"-'' ^ -N^u^'u'' + gf. (6) 
The familiar source-free inhomogeneous Maxwell Equations in physical spacetime remain 

WbHll = 0. (7) 

The more familiar form of the constitutive relations 

P[>c] = /ii^[t>c], (8) 

have been replaced in Eq. (|4]) by an equivalent single equation using Gordon's metric. Just as with a real optical 
metric (see [ll, 0, 0, Q ) , the source-free inhomogeneous Maxwell equation fT]) can be rewritten as 

= 0. (9) 

The covariant derivative in Eq. ^ is with respect to the complex optical metric, and except for the reciprocal of the 
impedance, = ^JTf\x^ would be the same as Maxwell's vacuum inhomogeneous equations in the optical spacetime 
but without polarizable materials, i.e., the components of F^ were obtained in Eq. ([4]) by simply raising indices on 
^ofc" using the optical metric. Because e and /i are ordinarily wavelength dependent, the values of N and Z at a 
spacetime point depend on the particular geometrical optics wave being considered. 

To proceed further with the geometrical optics approximation we must make assumptions about the size of the 
imaginary part, k, of the index of refraction in Eq. ([1]) at the frequency of interest. We will consider two types of cases; 
for the first type the imaginary part k is not small compared to the real part n (at the wavelengths of interest) and 
for the second type it is, i.e., k << n. The first type includes the absorption of low frequency waves in a conductor 
as well as the absorption of microwaves by water. The second includes a case of interest to us, the extinction of light 
waves traveling in a dilute gas. In the first case the eikonal S in Eq. Q has an imaginary part which is not negligible 
compared to the real part (see Sec.|Tl| and in the second type, S can be taken as real (see Sec.lIIip. In Sec.|n]we include 
two complex eikonal examples and in Sec. lIIII we give one real example which we then use to evaluate distance-redshift 
in standard cosmologies when absorption is present. Because much of the algebra of the second type is included in 
the first we start with a complex S. We give some concluding remarks in Sec. lIVI 



II. A COMPLEX EIKONAL 



In this section we develop the geometrical optics approximation for waves traveling in a medium where absorption 
on the scale of a wavelength cannot be neglected. Such significant absorption requires the use of a complex eikonal. 
By inserting Eq. ^ into Maxwell's Eqs. ^ and © we find that Aab is of the form 



Aab = -'2'k[aEb], 



(10) 



3 



where ka = daS is a complex null vector of the optical metric {k°'ka = 0) satisfying a complex geodesic like equation 
k = in the optical spacetime. In general, the invariant derivative '''is defined by 

'■' = ~k'Wb, (11) 

and rather than being a directional derivative as it is when the metric is real, it becomes a complex partial differential 
operator. 

Because we are interested in "homogeneous" waves, i.e., those for which the surfaces of constant phase and constant 
amplitude coincide we require that the spatial part of fc", i.e., the part of fc" which is orthogonal to u'^, be 
proportional to a real unit spacelike vector fc", k^'ki, = 1. This direction defines the propagation direction for the wave 
as seen by the optical fluid. As a consequence we can write 

r = -{S,bu'')N{Nu'' + 

ka = -{S.bu''){Ua+Nka). (12) 

The local period T and decay time of the wave are related to changes in the real and imaginary parts of the eikonal 
as seen by an observer moving with the fluid, i.e., by 

-(5.;,"'') = Ao(^^-*Jr), (13) 

which are in turn related to the local wavelength A and absorption coefficient a through the complex index of refraction 
Nhy 

27r a / 27r 2 \ 

-+^2-^(^-^^)- (^4) 

Maxwell's equations further restrict the electric field amplitude £a in Eq. (fTO]) by Eak"" — 0, but leave the remaining 
freedom of definition £a £a + f{x)ka (here f{x) is an arbitrary complex function). The first order (polarization 
dependent) correction to geometrical optics is given by 

Sab = 2(f[,,,] -fc[,I?b]), (15) 

with a remaining freedom Va — > "Da + g{x)ka- Furthermore the propagation equation for the electric field amplitude 
f " is 

t + £"9 + f > = y (Vb£'' + kbV'' + 24>,b£^), (16) 

where 2(/> is the natural logarithm of the reciprocal of the impedance, i.e., 20 = log -^/e//! , and 9 is defined as the 
divergence of fc". It is a generalization of the expansion rate of the "null rays" defined by the complex vector field fc°, 
i.e., 



and for a real metric is conventionally interpreted as the fractional rate of change of the observer independent 
cross-sectional area A of a small beam of neighboring rays. In what follows we choose a gauge where fait" = which 
makes £a spacelike and transverse to the wave's propagation direction, i.e.,£ofc° = 0. By contracting Eq. (fTH]) with 
£a we arrive at the propagation equation for the amplitude of plane polarized waves 

(£>) + 2(£,f'')(0 + 0) = 0, (18) 

which can be simplified to read 

\£a£'')A./7/^] = 0. (19) 



The time averaged 4-flux seen by an observer moving with the optical fluid is in general 
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S'' = —Re <j H'^'F.b 



(20) 



where {•} stands for 'the complex conjugate of. 

When F°-^ in Eq. @ is restricted to the lowest order geometrical optics approximation, Eq. ([2]), and is homogeneous, 
i.e., satisfies Eg. p^ . we have 



&|2t 



Re{7V} + k" 



(21) 



where 5/ is the imaginary part of the eikonal 5, and Re{A^} is the usual index of refraction, see Eq. ([T]). The coefficient 
of the fluid velocity in 5*° is the time average of the energy density ( x c) and the coefficient of k'^ is the time average 
of the magnitude of the Poynting vector, both measured by observers moving with the optical fluid. Equation (I2ip 
shows that energy in the single frequency geometrical optics wave is transferred in the k'^ direction with a speed of 
c/n by this wave. In the next two subsections we give two concrete examples where S is complex. 

A. Plane Waves in Minkowski Spacetime 



To make contact with familiar examples in classical electrodynamics, we start with a plane wave propagating in 
an optical fluid which is at rest in flat spacetime. We assume e and /i have only a z-dependence, and study waves 
propagating along the z direction, starting aX z — —oo. To suppress reflections and to make the geometrical optics 
approximation valid, we assume that e anf fi vary slowly over a wavelength A, i.e.,e^zA <C 1, /1.2A <C 1. The physical 
metric is flat Minkowskian, the fluid's 4- velocity is = ^o, and the optical metric, Eq. ([5]), is 



[cdtf 

' N{zy 



+ dx^ + dy^ + dz"^. 



(22) 



From Eq. (I12p we flnd the complex wave vector 



with the complex eikonal 

S = Sb. + iSi 
Equations and reduce to 



k"- = 7V(7V, 0,0,1), 
ka = (-1,0,0,7V), 



-ct+ / n{z')dz' 



K{z')dz' 



(23) 



(24) 



27r a ,,/27r 

= N { — 

A 2 \cT 



(25) 



which give the wave a constant frequency v = 1/T = c/(27rAo), but a z dependent wavelength A = 27rAo/n(z) and a 
z dependent absorption coefficient a = 2k{z)/Xo. For this wave we have chosen the scale of S so that the geometrical 
optics expansion parameter Aq is the wave's rationalized wavelength in the absence of refractive material. This wave 
corresponds to a constant and uniform source (at z — —00) which resulted in = 00 in Eqs. (jl3[) and (jl4p . 

The expansion defined in Eq. (|17p vanishes for this plane wave, i.e., 9 = 0. If it is linearly polarized along the x 
direction the amplitude of the field is (0,5^^,0,0). The propagation equation (fT9| then simphfies to 



[f-(e//i)V4] = 



which implies 



£^ ( - 



1/4 



f[-ct+ Niz')dz'], 



(26) 



(27) 
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where the function / reflects the time dependence of the source amphtude £{t, z) at the source, i.e., at z = — oo. For 
a stable plane wave source, we simply put / = constant. 

To compute the energy flux we can use either the spatial part of the general result Eq. ([2T]l , or because the physical 
metric is flat and the fluid is at rest, use the familiar 3-D electric and magnetic flelds from Eqs. ^ and PU)) 

E = e^^/^«£^i, 

H = y/7fii k X E, (28) 
to evaluate the time averaged Poynting vector directly 

S = -^Re{ExH}, 

= — Re{V^}|f"|'e-2S^/^«k. (29) 
Stt 



The magnitude of S can be evaluated using Eq. (I27p as 



^ cos/3(z) 
cos/5(— oo) 

where S{—oo) is the flux at the source and where 

c„,,.M^.Ma^ (31, 

The phase (3 represents the angle by which H fleld lags behind the E field, and cos (3 is the familiar power factor in 
the language of circuit analysis. From the imaginary part of the eikonal we can now easily write down the relation 
between the absorption coefficient a and the classical optical depth t Q 

t{z) = 2- / K(z')dz' ^ I a{z')dz'. (32) 



B. Spherical Waves in a Static Spherically Symmetric Spacetime 

Because this case is similar to the above, we truncate the discussion and give mainly the results. For an isotropic 
monochromatic source at rest at the origin of a static spherically symmetric spacetime which is emitting radiation at 
a steady rate into an optical fluid, which is also at rest, and whose optical properties depend only on the distance 
from the origin, we have an optical metric of the form 



g2<E.(r) 

' N{r) 



;{cdt)^ + e^^<-'-Ur^ +r^de^ + 



sm 



(33) 



and a fluid at rest with respect to the non-rotating Killing flow, i.e., = e *(5g. For the radial null vectors in Eq. ([12 
we have, 



ka - (-l,Afe-*+*,0,0). 



(34) 



The complex eikonal is 



S 



-ct 



n{r )e 



+ i 



K,{r )e 



(35) 



Equations (fT3|) and p4)) simplify to 



A 2 I cT ' 



(36) 
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and give an r dependent frequency ly = 1/T = ce */(27rAo), an r dependent wavelength A = 27rAoe*/n(r) and an r 
dependent absorption coefficient a — 2K(r)e~*/Ao. The expansion parameter 9 of Eq. (|f 7p is 



For a time independent source Eq. ()f 9p now gives 

[r(e/M)i/45] =0, (38) 

where the polarization vector has been written in the form f — (0, 0,£/r, 0) and points in the eg direction. The flux 
measured by u"" = e~*^Q is found from Eq. (|2ip to be 

on 

cos/3(0) W ^''''^ 

where C is the total isotropic power radiated by the stationary point source in a narrow frequency range. The optical 
depth T changes from Eq. to 

r(r) = 2Ao"' / K(r')e-*+*d/, (40) 
a(r')e*dr'. (41) 



The spherical result differs from the plane wave result of Eq. (|5D|) by a decrease of the flux caused by the wave's 
expansion {A~^ cx r~^) and by a frequency shift in the wave as it moves through the changing gravity field. 



III. A REAL EIKONAL 



Waves traveling in spacetimes where the imaginary part of the index of refraction is much smaller than the real 
part must still satisfy the same set of Maxwell's equations ([3]) and ([9]) as before but for them the eikonal S in Eq. ([2]) 
can be taken as real. The complex index of refraction is caused by a complex permittivity e = en + iej and/or a 
complex permeability = + i^j whose imaginary parts are much smaller than their real parts. Consequently we 
write N as 

N = n + iK = n{l + iXaTi + 0{Xo^)) , (42) 

where n = y'enfiR and the constant parameter Aq is the same parameter used to keep track of the various orders in 
the geometrical optics expansion. The absorption coefficient a of Eq. (|14p is related to k by a = 2kAo/A. The optical 
metric components of Eqs. ([5]) and ([6]) are 

gab = -TLuatifc + .g_Lafc = 5afc - 2 — -^Malife + O(Ao^), (43) 

-ab ^ -N^u'^u'' + gf = g'''' + '2^Kn'^u''u'' + 0{Xo'^), (44) 
where the ©(Ao") term of the optical metric, gab, and its inverse, g'^^, are real 

gab = ^UaUf, + g^ab, (45) 

g"'' = -n'^u'^u'' + gf. (46) 



When the geometrical optics expansion of Eq. ^ is inserted into Eqs. ([3]) and ^ the 0(Ao ^) terms result in Eq. ([T 
again, but now with ka the gradient of the real eikonal S and null with respect to the ©(Aq") optical metric gab, i-e., 

fc'fc, = 0, (47) 
= ^ = fc^V^r = 0. (48) 
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This gives real geodesies x°'{£) with real tangents dx°'/d£ = = g°'''ki,, and makes the real metric gab the important 
geometric quantity rather than the complex gab- The ' ' ' derivative is now the familiar derivative with respect to 
an affine parameter £ along null geodesies. Equations (fT2|l are still valid except TV is replaced by its real part n and 
complex fc° by the real fc". Equations (fT3|) and (fT4|) for the frequency and wavelength are replaced by 

- {S^bu") = Ao (^) , (49) 



cT 



and 



27r /27r\ , , 

The polarization vector £a in Eq. (fTO|) is now constrained to be orthogonal to fc" and can again be chosen orthogonal 
to the fluid Ua- The 0{Xo^) correction terms Bat are still given by Eq. the expansion 6 in Eq. ([T7]) is computed 
using k"", and the covariant derivatives are all with respect to the real gab Christoffel connection. The only equation 
that contains a new term, the extinction term, is the propagation equation for the amplitude f that replaces Eq. ([T| 



t + re + E"^ + ■Kn'^{u'^kdfE'' = y(V6^'' + kbV^ + 20,f,£^). (51) 

The phase is now computed similarly as for Eq. (|16p but now using only the 0{Xq^) terms of the impedance 
24> = log \/ €r/ IJ.R- The integral of Eq. ((5T|) which replaces Eq. p9|) now contains an affine parameter integral 



log [(£a^")AVefl/MiiJ = -2y Kn2(M^fcrf)2d£. (52) 
When the time averaged Poynting vector, Eq. (|20p . is evaluated Eq. (PTjl is replaced by 



S'^ = —{£a£''){S^bv'f (VWm^) + fc"J ■ (53) 

Comparing Eqs. pT|) with (|53p the effects of absorption can easily be seen to have shifted from the eikonal where 
it belongs if absorption is significant over wavelength scales, to the slowly changing amplitude £a where it belongs if 
extinction is significant only over many wavelengths. 

In the next section we give an example of absorption with a real eikonal from observational cosmology. We evaluate 
luminosity-distance in a Friedman-Lemaitre- Robertson- Walker (FLRW) universe when extinction occurs. 

A. Robertson- Walker Spacetime 

As an example of weak absorption we apply our complex extension of Gordon's optical theory to a Robertson- 
Walker (RW) universe filled with a time dependent index of refraction N{t) = n(i) -I- iK{t) [see Eq. (|42p ] and obtain 
for the real optical metric, Eq. (HS)) . 



= + + r\de^ + sin^ 1 (54) 



i{tY ^ ' i 1 - fc 

The familiar curvature parameter k — (1,0,-1) distinguishes respectively between spatially closed, flat, and open 
models. The radial outgoing null geodesies can be solved immediately to give null vectors 



and the related real spherically symmetric eikonal centered at the emission point, t = t^, r = 0, 

^(.,.) = 5.^i?o(-£^+sinn-H), (57) 
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where 

{sin[r] k = +1, 
r k = 0, (58) 

sinh[r] fc = — 1. 

The constant Rq is the current radius of the universe [from Eq. ([51)1 ] and has been introduced so that the geometrical 
optics expansion parameter Xq corresponds to the rationahzed wavelength of the wave when it reaches an observer at 
to,ro from an emitting source at te,r = 0. From Eqs. ([^9]) and ([50)) we have the wavelength and frequency redshifts 

l + Zn^^ = ^. (59) 



which are thus related by 



(l + z„) = ^(l + z). (60) 
n 



The conventional RW redshift {1 + z) = Rq/R is valid for both wavelength and frequency when the refracting and 
absorbing material is absent, i.e., when N = 1. We see that even with refraction and absorption the wavelength 
redshift remains as in RW cosmology. The frequency redshift, however, is affected by the real part of the index of 
refraction, n, but not by the imaginary part, n, i.e., not by extinction. 

To evaluate the apparent brightness of a source we need to evaluate the magnitude of the spatial part of the 
Poynting vector, Eq. (|53p . which requires that we know the area A in Eq. (|52p . For a spherical wave emanating from 
the comoving origin, A oc {rRY , which is confirmed by evaluating the expansion parameter 9 using Eq. (|55p . From 
the transport Eq. ([5^ we have 



{Sj-).J7^,^^j^, (61) 



where the optical depth T[t) from t to io is 



= 2c/ , (62) 



c/"^..'. (63) 



When looking back in time from the observer, Eq. (|53p gives the magnitude of the spatial part of flux 

°^ WmWmm?- ^^^^ 

The apparent luminosity L in a narrow band is just the value of S{t) at the observer, t = to- If the absorption is 
frequency independent in the observed frequency range, and the constant of proportionality is evaluated at the source 
at r{te) — the luminosity becomes 



'^TrRlrl {1 + Znf 



.-o2^2 n \ . (65) 



where C is the absolute luminosity of an assumed isotropically radiating source in that frequency range. The luminosity 
distance is easily read from this and differs from RW cosmology, i.e., 

dL = {l + z„)Roroe-'l^, (66) 

however, the apparent size distance_p| remains — r^Re- These distances obviously violate the classical reciprocity 
relation = (1 + zYdA (see e.g. Q) which is valid for non-refractive non-absorptive optics. For a discussion of the 
impact of violating the reciprocity relation on cosmology see Bassett & Kunz [10| . 
Using the dynamics of the FLRW cosmologies, we find 

c e^(^")/2 . r^'-'> dz' 



dL{Zn) = {1 + Zn)jJ y=-fimn^/\^l^\ (67) 

Hq Jo n{z')h{z') 
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where we have used Eq. (|60|) to express distance- redshift as a function of the frequency redshift z„, where Hq is the 
Hubble constant, and where 



h{z) = y/flA + W + z)2 + n„^il + z)^ + + z)^ . (68) 

The density parameters f^A, are standard and represent current relative amounts of non- interacting gravity 

sources: vacuum, pressureless matter, and radiation energies. They are related to the curvature parameter fi^ by 



When the optical depth is written as a function of the frequency redshift z„ it becomes 



^k^ —f7TB2 = l- {^A + ^rr. + ^r)- (69) 

-tlnHr, 



Hq Jo (1 + z')n[z')h[z') 

The distance-redshift of Eq. (|67p can be compared with a similar result but without absorption given in Q . In that 
paper we have shown that a cosmological model with a refraction index n{z) = 1 + az^ + bz^, but with no absorption, 
can fit the currently available type la supernovae data quite well. A source for such refraction remains elusive. For 
illustrative purposes, we now use the distance-redshift of Eq. (|67p with a constant absorption coefficient a which does 
not change with frequency (i.e., is grey) but without refraction to fit this same supernovae data, i.e., we take 

r^A = 0, n{z) = 1, a{z) = const. (71) 

The Hubble constant we use is Hq — 65 km • s~^ • Mpc^^. Since we are concerned with the matter dominated era, we 
exclude radiation {fir = 0). The fimction h{z) in Eqs. (|57|) and ([70)1 simplifies to 

h{z) = {l + z) ^/^ + ^lmZ, (72) 

and our model has just two parameters (rJm, en). We compare the distance modulus versus redshift, /i(z), of the 
concordance model (J7a = 0.7, ftm = 0.3) with five JIa = extinction models (three open, one flat and one closed). 
The result is shown in Fig.[Tl where the distance-modulus /x, defined by 

^ = ^l°Sl^+2^' 

is compared to the supernova data [ll|, [l^, [3i Hj We use the 178 supernova from the gold sample [l^ with 
redshifts greater than cz = 7000 km/s. The critical redshift region is between 0.2 < z < 1.2, where most of the 
supernova data is concentrated. In Fig.iU the solid green curve represents the flat concordance model dominated 
by dark energy. The dotted red curve is an open universe with baryonic only matter, fl^ = 0, fim = 0.05, n = 1 
and a{z) — 7 x 10~^ Mpc~^. The short dashed red curve is an open model containing only dark matter, — 0, 
fini = 0.3, n = 1 and a{z) — 1.3 x 10"'' Mpc"'^. The solid black curve is our least model (see Fig. [5]) with fl\ — 0, 
flni — 0.73, a{z) — 2.2 x 10^** Mpc^'^. The longer dashed green curve represents a flat universe dominated solely by 
matter, with — 0, — 1-0, a{z) = 2.6 x 10"'' Mpc"'^. The longest dashed blue curve is a closed model, where 
Ha = 0, = 1.3, a = 3.2 X lO"" Mpc'^ The black dashed curve is the now disfavored matter dominated J^a = 0, 
^^m — 0.3 model (without extinction). In the inset we show Afi versus z (relative to the disfavored matter dominated 
case) for each model. To produce roughly the same amount of change in the distance-modulus a larger absorption 
coefficient a is needed for a larger {a is positively correlated to flni)- As the reader can easily see in Fig. [1] the 
effects of a suitable value of the simplest absorption coefficient a{z) = const can simulate the accelerating effects of a 
cosmological constant. 

In Fig. 121 we show the confidence contours for our model parameters, a versus flm- The best fitting parameters 
are flm = 0.73, a = 2.2 x 10~'''Mpc^^, with Xmin = 1-04 (per degree of freedom). The innermost contour encloses 
the 68.3% confidence region, and the next one encloses the 95.4% confidence region, the outermost one encloses the 
99.73% confidence region. For each fixed fl^ selected in Fig.|l] a is the value that gives the least x^. 

If we extract a density factor p from the absorption coefficient, i.e., a = a/ p we obtain an opacity in e.g., cm^ • g""'. 
The density p with units of g • cm~'^ is the density of the relevant species causing the absorption. A competitive 
absorption model should properly account for both the cosmic expansion and the physical/chemical evolution of 
the inter-galactic medium [17|. Here we are content with an order of magnitude estimate noting that an opacity 
a = 10^ cm^-g^^ as proposed for the carbon needle model in Aguirre [l^, requires a density p the order of 10""^^ g-cm~^ 
to produce the absorption needed for the above fitting. This density is only a factor of ~ 10~^ of the current critical 
mass density pc ~ 8 x 10~^°g • cm~^. For fine-tuned dust absorption models, see e.g., p^. [l9l [20l. [2ll. [22j . 



10 




0.0 0.5 1.0 1.5 2.0 



z 

FIG. 1: Distance modulus fi versus redshift z. The two red and one solid black curves are open models, the two green curves 
are flat models, and the blue curve is a closed model, (looking downward at redshift z — 1.5) Dotted red curve: JIa = 0, 
= 0.05, a{z) = 7 X 10"^ Mpc"^ Short dashed red curve : Qa = 0, = 0.3, a{z) = 1.3 x 10"" Mpc"\ Solid black curve 
(our best fit): Oa = 0, f^m = 0.73, a{z) = 2.2 x 10"" Mpc"\ Solid green curve (concordance model): JIa — 0.7, Qm ~ 0.3, 
Of = 0. Longer dashed green curve: Q,a = 0, Q,m = 1.0, a{z) — 2.6 x 10"" Mpc"^. Longest dashed blue curve: fim = 1.3, 
a{z) = 3.2 x 10"" Mpc"\ Black dashed curve: 57a = 0, flm = 0.3, a = 0. Inset: Afi versus z curve for each model, the fiducial 
model (black dashed curve): 57a = 0, fim ~ 0.3, a = 0. 
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FIG. 2: a versus Qm- The best fitting parameters are fim ~ 0.73, a = 2.2 x 10~"Mpc"^, with Xmin = 1-04 (per degree of 
freedom). The innermost contour encloses the 68.3% confident region, the next one encloses the 95.4% confidence region, and 
the outermost one encloses the 99.73% confidence region. 
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IV. DISCUSSION 

We have demonstrated how the 4-D optical metric of Gordon [l| [see Eq. ([5])] can be extended and used even in 
cases where absorption is present. We looked at both 'strong' and 'weak' absorption, the distinction being whether 
absorption is significant on wavelength scales or only over a multitude of wavelengths. The two cases are distinguished 
respectively by complex and real eikonals. For the complex eikonal case the optical metric must remain complex (see 
Sec.|lT|), however, for the weak absorption case the real part of the optical metric (essentially the same as Gordon's 
original proposal) remains as the significant geometrical structure (see Sec. lIII|) . The two cases differ on how absorption 
appears in the geometrical optics field. In the 'strong' absorption case the imaginary part of the eikonal reduces the 
wave's intensity [see Eq. (|2T|) ] but in the 'weak' case the amplitude's reduction [see Eq. ([52]) ] is responsible for the 
intensity decrease. 

A geometrical optics wave is like a single frequency wave even though the wave's frequency changes from spacetime 
point to spacetime point. To superimpose multiple frequencies is straightforward, however, to superimpose optical 
metrics, real or complex, makes no sense. Consequently a useful single optical metric only exists when the optical 
properties are insensitive to the superimposed frequencies. Such frequency independence approximations are often 
designated as 'grey' in astrophysical applications. 

The two examples we gave in Section [TT] can be used to study the impact of refraction and /or absorption on light 
propagation in stellar atmospheres. The classical radiation transport equation (see e.g. [23l, [24|) is derived assuming 
that light follows null geodesies in curved spacetimes. With the presence of light refraction, both the direction and 
speed of light change. This would require the radiative transfer equation to be written using the optical metric instead 
of the physical metric. The optical metric might also be of some use in hydrodynamical simulations of stellar interiors. 
For an example, the slowing down of light will reduce the efficiency of energy transport outward via the radiation 
field. Our first example (flat spacetime) can be used to study the 1-D case; a comparison of results of a numerically 
solved radiation transport equation with/without a refraction index n(z) would be interesting. Similarly, our second 
example (curved spherically symmetrical spacetime) can be used to model the atmospheres of neutron stars. Further 
non-spherically symmetric examples could be useful in radiation transport calculations in accretion disks of black 
holes. 

Our last example, the real eikonal case, allowed us to give luminosity distance redshift for observations in standard 
cosmologies where both refraction and absorption were present. As an example of the usefulness of this theory we 
went on to fit the gold sample of type la supernovae to a Hubble curve corrected for grey extinction. In Section IIIII 
we showed that it is possible to explain the current supernova observations via a simple absorption model instead 
of requiring the existence of dark energy. Our best fit was an open = 0.73 model with constant absorption 
a — 2.2 X 10^'*Mpc~^, see Fig. [2] for confidence contours. More realistic z dependent models for a are in order. 
Since the flat concordance model is supported by other observations, e.g., cosmic microwave background data and 
baryonic acoustic oscillations, an absorption theory cannot be on firm ground unless it provides an explanation for 
these additional observations. We leave these and other applications to future efforts. 
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